Geometrodynamics in a spherically symmetric, static crossflow of null dust 
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The spherically symmetric, static spacetime generated by a crossflow of non-interacting radiation 
streams, treated in the geometrical optics limit (null dust) is equivalent to an anisotropic fluid 
forming a radiation atmosphere of a star. This reference fluid provides a preferred / internal time, 
which is employed as a canonical coordinate. Among the advantages we encounter a new Hamiltonian 
constraint, which becomes linear in the momentum conjugate to the internal time (therefore yielding 
a functional Schrodinger equation after quantization) , and a strongly commuting algebra of the new 
constraints. 



O 
O 

(N 
-t— > 
O 

O 

(N 



> 



in ■ 

o : 
\o . 
o ■ 

cr 

i ' 



X 



I. INTRODUCTION 

Beside the covariant description of general relativity, a 
Hamiltonian formalism of gravity based on the existence 
of foliations for globally hyperbolic space-times was de- 
veloped by Arnowitt, Deser andMisner (hereafter ADM) 
P|. In this approach the canonical coordinates are the 
components of the induced metric on the 3-leaves of the 
foliation, while the canonical momenta are related in a 
simple way to the extrinsic curvature of these spatial hy- 
persurfaces. Gravitational evolution is therefore quoted 
as geometrodynamics. The freedom to perform coordi- 
nate transformations on the leaves of the foliation leads 
to the diffeomorphism (or momentum) constraints. The 
true dynamics is encompassed in the so-called Hamilto- 
nian constraint. These four constraints (per point) form 
a Dirac algebra 2], which is not a true algebra in the 
mathematical sense, as its closure is obstructed by the 
appearance of the induced metric in the Poisson brackets 
of the Hamiltonian constraints. 

The problem of time in canonical gravity was reviewed 
by Isham in Ref. |3J- Approaches for introducing the 
concept of time are of three types: time is either identi- 
fied before or after quantization, or in certain approaches 
time plays no fundamental role at all. In what follows, 
we are interested in identifying time at the classical level. 
Time is not preselected by any Hamiltonian description of 
gravity, there are infinitely many ways to choose the time 
(many- fingered time formalism) Despite this ambi- 
guity, in certain cases it is possible to select a preferred 
time function, by cither imposing coordinate conditions 
[f| or by filling space-time with an adequate reference 
fluid [|| , and letting gravity to evolve in the (proper) 
time of the chosen reference fluid. 

In certain cases new canonical variables can be intro- 
duced, providing new constraints for gravity @, @, [Toj - 
Among the advantages we count that the Dirac algebra 
transforms into a true algebra and the quantization of the 
Hamiltonian constraint, usually leading to the Wheeler- 
de Witt equation (which has no linear space of solutions), 
rather gives a Schrodinger equation. This program has 
been particularly successful for incoherent dust, as pre- 
sented by Brown and Kuchaf in Ref. [|J. 

A similar formalism p"l| was by Bicak and Kuchaf 



applied for null dust, the geometrical optics approxi- 
mation to non-gravitational radiation. Null dust how- 
ever provides no natural time function, basically because, 
unlike the congruence of the incoherent dust particles, 
null world lines have no natural parametrization. While 
for ordinary dust the Hamiltonian- and supermomentum 
constraints depend on four pairs of canonical variables 
associated with the proper time and the comoving coor- 
dinate frame of the dust, the constraint equations for null 
dust contain only three pairs of comoving coordinates. 

The quantum theory of gravitational collapse can be 
modelled in the most simple sphericall y sy mmetric case 
by a collapsing thin shell of null dust |12|. [l3| . A sec- 
ond null dust shell can be introduced in the model in 
order to test the quantum behaviour of the geometry 
induced by the first shell. Motivated by certain prob- 
lems in the above scenario, the canonical formalism in 
the presence of a null dust has been recently extended to 
the case of two cross-flowing, non-interacting null dust 
streams in a spherically symmetric space-time by Bicak 
and Hajicek |14| . This formalism combines ingredients of 
the canonical formalisms developed for null dust with 
elements of the geometrodynamics of the Schwarzschild 
space-time [l^, developed by Kuchaf. The lack of a 
time-standard for a single null dust however deprived 
the canonical formalism of the cross-streaming null dust 
from a time-standard as well. This is because the start- 
ing point of the canonical description |l4j is simply the 
sum of the spherically reduced Einstein-Hilbert action for 
gravity and two pieces of the null dust action, also re- 
duced by spherically symmetry. The null dust variables 
are therefore doubled, without any of them becoming an 
internal time. The basic assumption of Ref. is that 
the cross-flowing null dust streams interact only gravita- 
tionally, therefore the energy-momentum tensors of the 
components are conserved separately. The analysis of 
the equations of motion provides two pairs of integrals of 
motion (per point), one pair for each null dust compo- 
nent. Unfortunately the Hamiltonian density could not 
be explicitly expressed in terms of these quantities, ex- 
cept in the case when one of the null dust components is 
switched off. In this case the action can be transformed 
such that the matter part of the Liouville form contains 
the integrals of motion associated to the null dust com- 
ponent in question. 
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The formalism derived in Ref. (lif is valid for cer- 
tain known spherically symmetric space-times, for exam- 
ple the Vaidva space-time, describing the one-component 
null dust |l6|. and the static space-time found in Ref. 
[l7| by one of the present authors. The latter space-time 
represents the geometry in the presence of a static cross- 
flow of non-interacting null dust streams. Although it is 
asymptotically non-flat and it has a central naked singu- 
larity, it can be conveniently interpreted as the radiation 
atmosphere of a star. A second interpretation presented 
in Ref. 01 is °f a 2-dimensional dilatonic model, in the 
presence of a pair of 2-dimcnsional scalar fields. While 
the dilaton is the square of the radial coordinate, the 
scalar fields are related to the energy densities of the 
null dust streams. The third interpretation, based on 
previous work of Letelier |l8j | , is of an anisotropic fluid, 
with radial pressure equal to its energy density and no 
tangential pressures. The static solution nas a n0 ~ 
mogenous counterpart pjj, which can be interpreted as 
a Kantowski-Sachs type cosmology. These two space- 
times obey a unicity theorem, as they are the only spher- 
ically symmetric solutions of the Einstein equation in the 
presence of a cross-flow of null dust streams with an ad- 
ditional (fourth) Killing vector 19]. Interestingly, for 
null dust streams with negative energy density, worm- 
hole space-times emerge [2(|, |2l| . 

The anisotropic fluid interpretation of the static space- 
time the cross-flow of null dust streams with positive 
energy densities is particularly important for our pur- 
poses. The physical model of the anisotropic fluid has 
a preferred time, which is the time elapsed in the rest 
frame of the fluid. This suggests that in contrast with 
the single null dust model, for the two component null 
dust an internal time formalism can be constructed. In 
this paper we will explicitly construct the matter action 
for the static configuration of non-interacting null dust 
streams in terms of suitable variables, containing the in- 
ternal time singled out uniquely by the cross-flow of null 
dust. 

In Sec. II we summarize the basic ingredients necessary 
for the purposes of the present work. We present: 

(A) the canonical formalism of ordinary incoherent dust 

|6j, with special emphasize on how the proper 
choice of the internal time allows us to introduce 
a set of new constraints for gravity, such that the 
new super-Hamiltonian constraint becomes linear 
in the canonical momentum conjugate to the inter- 
nal time; 

(B) the geometrodynamics of the spherically symmetric 

static vacuum |15| . with special emphasize on the 
introduction of geometrically motivated canonical 
variables (including the Schwarzschild mass) in the 
gravitational sector; 

(C) the spherically symmetric, static space-time with 

crossflowing null dust streams |17| and 



(D) the anisotropic fluid interpretation of the cross-flow 
of non- interacting null dust streams |18| , which pro- 
vides the internal time for the two component null 
dust. 

In Sec. Ill we introduce an action functional of three 
scalar fields characterizing the static cross-flow of null 
dust minimally coupled to gravity. We show that varia- 
tion with respect to the metric together with the equa- 
tions of motion reproduces the energy-momentum tensor 
of two non- interacting radiation streams. Two pairs of 
conservation equations for the rest mass currents and the 
momentum currents also emerge. 

In Sec. IV we derive the contribution of the two null 
dust streams to the super-Hamiltonian and diffcomor- 
phism constraints. Then we fulfill the program of re- 
placing the total super-Hamiltonian and diffeomorphism 
constraints by an equivalent set, in which both momenta 
conjugate to the temporal and radial canonical variables 
appear linearly. We also prove that the new constraints 
form an Abelian algebra. 

Sec V. contains a discussion of the falloff conditions 
the gravitational variables, the lapse and the shift should 
obey. 

In Sec. VI we compare our findings with the results 
presented in Ref. [14| and we show that similar tech- 
niques can be employed in the more generic context of 
Ref. as well. We also underline the connections be- 
tween our canonical variables and those employed in Ref. 
fl4|. specified for the static case. 

Finally in Sec. VII we summarize our results. 



II. PRELIMINARIES 

In this section we present a more technical summary 
of the results of Refs. 

H, E3> E3 and 

needed later 

on in the paper. 



A. Geometrodynamics of space-times with 
ordinary dust 

The space-time action of ordinary dust was con- 
structed by Brown and Kuchaf from eight scalar fields 
Z k ,W k ,T,M (k = 1,2,3) minimally coupled to the 
space-time metric ^g a b- 



S 



D 



T, Z k ,M,W k ;^ g a 



= -\J d^x^i^MiUaU + 1) , (1) 

The four-velocity U a is expressed as the Pfaff form of 
seven scalar fields. 



U a 



-T a + W k Z\ a . 



(2) 
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The equations of motion are 
. SS 1 







SM 
SS D 
SW k 
SS D 
ST 
SS D 
SZ k 



= -V- {4) gMz k , a u a , 



= -(V^gMU a ), a 



= -{V- (4) gMW k u a ), a . 



(3) 
(4) 
(5) 
(0) 



According to Eq. (JIJ the three vector fields Z k are con- 
stant along the flow lines of U a (they can be interpreted 
as comoving coordinates for the dust.) Eq. © shows 
that the four-velocity U a is a unit time-like vector field. 
Eq. lO allows us to interpret M as the rest mass den- 
sity of the dust and it represents mass conservation. Eq. 

can be interpreted as the momentum conservation 
law. From Eqs. (J5J), (J3J) and J3J) it is straightforward to 
deduce that T is the proper time along the dust world 
lines, measured between a fiducial hypersurface T = 
and an arbitrary hypersurface with constant T . The dust 
energy-momentum tensor T ab can be found from the vari- 
ation of the action Q with respect to ^ 4 'g a b- From the 
conservation of T a j, and M it follows that the dust par- 
ticles evolve along geodesies. 

The Legendre transformed action is 



S D [T, Z k ,P,P k ,g ab ,N,N a ] 



= I dt I d 6 x{PT + P k Z k - NH U - N a Ha 



where g ab denotes the induced metric on the leaves, N 
and JV a are the lapse function and shift vectors, respec- 
tively, and the momenta P and P k are conjugate to T 
and Z k . (The original variables W k were expressed in 
terms of P and P k .) The constraints are 



H 



D 



1 P 2 



2 Mg 1 / 2 



1 Mg 1 ' 2 
2~P 2 ~ 



P 2 



„ct6 



H D HP) 



P k Z h , 



(8) 
(9) 



The dependence of the Hamiltonian constraint on the 
variable M is spurious. This can be shown as follows. 
By varying the action with respect to M we obtain an 
algebraic expression from which M can be given in terms 
of the other variables. Substituting this into the Hamil- 
tonian constraint gives 



H 1 



P 2 



g- b HPH° 



(10) 



so the mass multiplier M is eliminated from the action. 

By employing that the total (gravitational + dust) con- 
straints have to vanish, e.g. = — H± and HP = —Hf? 
on the constraint hypersurface, and solving the con- 
straints I|1U|) . @ with respect to the momenta, we can 
replace the old constraints by an equivalent set. The new 



super-Hamiltonian constraint can be cast into the form 
ff T = P + h[g ab ,p ab ] = , (11) 
-jijHpHf , (12) 



h 



.9' 



where p ab are the momenta conjugate to g a b- Similarly 
the new supermomentum constraint is: 

P k + h k [T,Z k ,g ab: p ab } =0 , 



H 



h k = Z a k H^ - hT a Z a k 



(13) 
(14) 

The quantization of the linearized constraint l|llfl gives a 
Schrodinger equation [6|. 



B. Geometrodynamics of spherically symmetric 
static vacuum 

After the preliminary studies on the canonical for- 
malism of the spherically symmetric space-times |23j . 
a comprehensive analysis of Hamiltonian dynamics for 
Schwarzschild black holes was given by Kuchaf [lq. In 
this section we summarize those results of his work which 
are relevant for our purposes. 

The space-time was foliated by spherically symmetric 
leaves St which were labelled by the parameter time t. 
The induced metric on these 3-leaves can be character- 
ized by two metric functions A and R, 



g ab dx a dx b = A 2 (i, r)dr 2 + R 2 (t, r)dQ 2 



(15) 



where r is a space-like coordinate and dQ 2 is the line 
element on the unit sphere. Under coordinate transfor- 
mations R behaves as a scalar and A as a scalar density. 
In the ADM decomposition of the spherically symmetric 
geometry, the shift vector has a non- vanishing component 
only in the radial direction, denoted with N r , which to- 
gether with the lapse function N depend solely on the 
variables t and r. 

The metric functions R and A are chosen as canonical 
coordinates and their momenta, as derived in ITF 



are 



Pa = —N~ 1 R(R — R'N r ) 
P R = -N- 1 \a(R - R'N r ) 



R(A - (AN r 



(16) 



The vacuum action for the spherically symmetric geom- 
etry can be written as 

Sx[g ab ,N, N a ] 

dt I d 3 x{AP A + RP R - NHf - N r H?) , (17) 

with super-Hamiltonian and supermomentum con- 
straints 

H±[R, A, P A , Pr] - ^PrPa + ^ AP 2 

+ LrR»-§R<A> + ±R> 2 - 1 -A, (18) 
H G [R, A, P A , Pr] = PrR' - P A A . (19) 
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There exists a canonical transformation, through 
which the only dynamical characteristics of the 
Schwarzschild space-time, the Schwarzschild mass M 
turns into a canonical variable. The new set of variables 
is (M, R; Pm, Pr), where M(t, r) is expressed in terms of 
the old variables (A, R; Pa, Pr) through the formula of 
the Schwarzschild mass derived by Kuchaf : 



m = \r~ x Pa - \^ 2 RR' + \r ■ (20) 



The remaining part of the canonical transformation is: 
2Af 



Pm = AP A 
R = R , 

Pr = Pr - 



1 - 



R 



R- 



2^ lAPA - 2 



1 / 2M \ 1 , , 



2M\ 
~R~ 



(AP A )' RR! - (AP A ) (RR')' 



(21) 



The second advantage of the new set of canonical vari- 
ables is that the momentum Pm is the gradient T' of 
the Schwarzschild time (cf. Eq. (80) in Ref. 0j). The 
gravitational constraints i|18fl and i|19|) . written in terms 
of the new canonical variables, become 



P^[M,R,P M) Pr] = - 1- — 



2M\ 1 M'R' 



R J A 
i _ 2M\ P M P R 



- + 



(22) 



R J A 

P r G [M, R, P M , P R ] = P R R' + P' M M , (23) 

where A rather than being a canonical variable, is only 
a shorthand notation for the following expression of the 
new canonical variables 



A 



2M 



M 



f/2 



2M\ 9 
1 - — P, 2 



R / 



M 



(24) 



We will also introduce the canonical variable M in 
the description of the gravitational sector of the cross- 
streaming null-dust space-time. 

We mention here the related result of Varadarajan [22| . 
who has derived a transformation from the usually em- 
ployed canonical variables (induced metric + extrinsic 
curvature), to a set of new canonical variables, which 
have the interpretation of Kruskal coordinates. This 
transformation is regular on the whole space-time, in- 
cluding the horizon. The constraints simplifies in such 
an extent, that those are equivalent to the vanishing of 
the canonical momenta. 



C. The spherically symmetric, static space-time 
with crossflowing null dust streams 

The static superposition of two non-interacting null 
dust streams propagating along the null congruence u a 
and v a is characterized by the energy-momentum tensor 



P 



2ND 
ab 



p{u a u b + v a v b ) 



with 



u a u a = v a v a = , u a v a £ 



(25) 



(26) 



The same time-independent energy density p was chosen 
for both null dust components in order to assure no net 
energy flow (static configuration). 

The spherically symmetric, static space-time contain- 
ing such a cross-flow of two non-interacting null dust 
streams has been presented in Ref. fl7| : 

L 1 

ds 2 = -2a-^— [dZ 2 - R 2 (L)dL 2 ] + R 2 (L)dn 2 , (27) 
R(L) 

where Z and L are the time and radial coordinates 
adapted to the symmetry and R is the following expres- 
sion of the radial coordinate: 



-R(L) = a(e L - 2L<t> B ) 
$ B = B + I e x ~dx 



(28) 



Here a is a positive constant and B is a parameter. 

The four-velocity null vectors of the null dust streams 
are then 



WZ, 



RWL, 



v n = WZ_ n - RWL_, 



with 



W = 



r 2 



R 



The energy density becomes 



p = (8itR 2 W 2 ) 1 



(29) 



(30) 



(31) 



The superposition of the in- and outgoing null dust 
streams can be interpreted as an anisotropic fluid. This 
indicates that there may be a possibility to use the 
same procedure as in the case of the incoherent dust 
to obtain an internal time for the canonical dynamics 
of cross-flowing (but otherwise non-interacting) null dust 
streams, minimally coupled to gravity. 



D. The anisotropic fluid interpretation of the 
cross-flow of non-interacting null dust streams 

Letelier has shown that the energy-momentum tensor 
of two null dust streams is equivalent with the energy- 
momentum tensor of a specific anisotropic fluid |18| . As 
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consequence, the source of the static, spherically symmet- 
ric space-times l|77|) can be interpreted as an anisotropic 
fluid with radial pressure equaling its energy density and 
no tangential pressures: 



We also note that 



Tab = p(U a U b + XaXb) ■ 



(32) 



Here \ a 1S the (normalized) radial direction and U a is 
the unit four-velocity of the fluid particles, obeying 

-u a u a = X aX a = i , u aX a = o . 

They are related to the null vectors by 

U<* = -±=(u a +v a ), x <* = -L(u a -v a ). (33) 

By employing Eqs. I|29ll . we can also express the vector 
fields U a and x a m the coordinate basis defined by Z and 
L: 



U a = s/2WZ, a , xa = V2RWL a 



(34) 



In the anisotropic fluid picture p represents both the en- 
ergy density and the pressure, while no tangential pres- 
sure components to the spheres of constant L are present. 
The fluid is isotropic only about a single point, the origin. 



III. ACTION PRINCIPLE FOR THE STATIC, 
SPHERICALLY SYMMETRIC CROSS-FLOW OF 
TWO NON-INTERACTING RADIATION 
STREAMS 

A generic spherically symmetric space-time, in coordi- 
nates (T, R, 9, ip) is characterized by two metric func- 
tions h and / as: 

ds 2 = -h(T, R)dT 2 + / -1 (T, R)dR 2 + R 2 dfl 2 . (35) 

Let us introduce two scalar fields Z (T) and L (R) , and 
the following advanced-type and retarded-type combina- 
tions of the 1-forms dZ and dL, which span the (T, R) 
sector: 



Ua = WZ^ + RWL^a 
V a = WZ M ~ RWL M . 



(36) 

with W given by Eq. (|30|l . Thus in this co-basis the 1- 
forms u a and v a do not have time-dependent components. 
They are entirely expressed in terms of the two scalar 
fields Z and L (as the coefficient functions W and R 
can be given in terms of L). Note that the expressions 
(|36[) are identical with (and in fact motivated by) Eqs. 
(|29[) . but this time the scalars Z and L are not related 
to any exact solution, and in consequence the 1-forms u a 
and v a are not necessarily null for the generic spherically 
symmetric metric (|35f) . They do have instead the same 
length: 



vTu a = v"v a 



dZ 

If 



fR 2 



dLY 

Irj 



(37) 



u v a 



-W 



dT 



(38) 



Let us define a dynamical system by the action: 

S 2ND [^g ablPl Z,L] = -~ f d 4 x^h)p(uau a +vav a ) , 

(39) 

where p (L) is a third scalar field. We do not know at 
this stage, what is the dynamical system described by 
the action (|39|l . 

Variation of the action with respect to the metric gives 
the energy-momentum tensor: 



rpab 



-2 SS 



2ND 



= ^g ab p(u c u c + v c v c ) + p(u a u b + v a v b ) , (40) 

while the variation with respect to the coordinates Z, L, 
and the parameter p give the Euler-Lagrange equations: 






= 



SS 



■2ND 



5Z 

SS 2ND 



SL 



= -2[V Z ^gpW(u a +V a )ia, 

= 2^-^)gpL{u a u a + v a v a ) 
-2R[^^)~gpW(u a ~v a )\ a , 



SS 



2ND 



iSn =-^V Z ^g(UaU a +V a V a ) 



(41) 

(42) 
(43) 



In Eq. (|42ll we have employed the relation dW/dL = 
W (2RL - dR/dL) /2R. 

Eq. lO together with Eq. implies that both u a 
and v a are null vectors. Then the energy-momentum 
tensor (|4Ti|> reduces to 



T 



ab 



p(u a u" + v a v") 



(44) 



characterizing a non-interacting cross-flow (in the null 
directions u a and v a ) of null dust streams with energy 
density p. 

One can define rest mass currents as in 



J a := V- {4) gpu a , JC a 
and momentum currents as 

Jl ■= WJ a , K a L := WK a 



gpv 



(45) 



(46) 



In term of these Eq. I|41|l is a continuity equation for the 
net flow of radiation: 

V a (J£+lC a L ) = 0. (47) 

As the vectors u a and v a are null, Eq. I|42() simplifies to 

V a {Jl-JCl) = 0, (48) 
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implying that both momentum currents are conserved 
individually: 



V a J£ = 



V Q /C£ = 



(49) 



as expected for non- interacting radiation fields. 

We have shown that the action defined by Eqs. i|36*|) 
and H39|) describes a cross-flow of non-interacting null 
dust streams in a static configuration with energy density 
p. As the vectors u a and v a are null, we can partially 
normalize them as u a v a = — 1. Also, from Eq. (|37[) we 

get dZ/dT = (fh) 1/2 RdL/dR. Then Eq. I® allows to 
express both metric functions as 



r 



2ae L R 



h = 



2ae L 



dL 
dR 

dZY 
dTl 



(50) 
(51) 



By inserting these into the generic spherically symmet- 
ric metric ()35fl. we obtain the metric form J2ZJ, however 
without the additional information 128fl and <|31|) . In or- 
der to recover these, we need the Einstein equations, de- 
rived from the sum of the Einstein-Hilbert action and 
the cross- flowing null dust action (l3*5|) . These are iden- 
tical to those presented in Ref. [l?J, thus lead to the 
solution summarized in Section II. C. 

At the end of this section we note that the equivalent 
action in the anisotropic fluid picture is 

S F [^g a \p,Z,L] = - 1 - J d 4 x^/^r g p(U a U a + XaX a ) , 

(52) 

with U a and Xa given by Eq. (|34|1 . Due to the equiv- 
alence of the two interpretations, all equations are the 
same, irrespective of they being derived from the cross- 
streaming null dust action 1)39(1 or from the anisotropic 
fluid action (JH21- 



IV. CANONICAL FORMALISM 

In this section we present the calculations yielding lin- 
earized constraints for the two-component null dust, sim- 
ilar to Eqs. and 1|13|) derived for ordinary dust. 



coordinates: p = p(t,r), Z — Z(t,r) and L = L(t,r). 
From Eq. ^ ^/Wg = N^/g . 

The (3+l)-split form of the Lagrangian density taken 
from the action Q39JI is 



L 



2ND 



gygpw 2 

N 

?2w2 



Z 2 -2N r ZZ'- 



N 2 - A 2 N r2 
A 2 " ' 



N 



A 2 



(54) 



The canonical momenta conjugate to the radiation vari- 
ables Z and L become 



Pz := 
Pl : 



dL 



2ND 



2a^gpW< 



dZ N 
dL 2ND ^ 2a^/gpR 2 W 2 



(Z - N r Z') , 

(L - N r L') , (55) 



or inverted with respect to the velocities we obtain 
N 



Z 
L 



2a^gpW 2 
N 



2a^/gpR 2 W 



Pz + N r Z' 



-Pl + N r L' 



(56) 



By inserting the velocities only in one factor of the 
velocity-squared terms of (|54(l we obtain the Lagrangian 
in the " already Hamiltonian" form 



-2ND 



ZP Z + LP,, - NH 



2ND 



r tj2ND 



N r H, 



(57) 



where the Hamiltonian and momentum constraints asso- 
ciated with the cross-flow of null dust streams are found 
to be 



H 



2ND 



1 



H: 



2ND 



2a^gpW 2 
2a^/gpW 2 
A 2 

Z'P Z + L'P L 



p2 

Fz+ R 2 



(Z 12 + R 2 L' 2 ) , (58) 
(59) 



Remarkably, the momentum constraint has the same 
form as the dust constraint ©. 



A. 3+1 decomposition of the two null dust 
Lagrangian 

The ADM decomposition of any spherically symmetric 
metric yields [l5|: 

ds 2 = -(N 2 -A 2 N r2 )dt 2 + 2A 2 N r drdt + A 2 dr 2 +R 2 dn 2 , 

where A and R are the metric functions from the induced 
line-element (|15H and (t, r) are generic coordinates or- 
thogonal to the (8, (p) sector . The variables p, Z, L char- 
acterizing the radiation cross-flow thus depend on both 



B. Introduction of new dust constraints 

If we vary the dust action (|54|l with respect to the 
comoving density p of the dust, we obtain 



XC2ND fifT^ND 

op op 
from which p can be expressed as 



_ n P 2 +Pf/R 2 

2a ^ w = A h> 2 + £ 2 L> 2 



(60) 



(61) 
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By substituting this result into the Hamiltonian con- 
straint we get 



H 



2ND 



(62) 

Since l|59l) implicates that the last two terms below the 

root appear in (H^ ND \ , we eliminate them from (1621) . 
The final form of the Hamiltonian constraint is 



H 



2ND 



We note that in the spherically symmetric case the mo- 
mentum constraint H 2ND can also be brought to a square 
root form. From Eq. I|15|) and l|63[l we have 



H 



2ND 



PzL'R 



P L Z' 
R 



(AH 



2ND\2 



(64) 

Eq. H63|) i s °f similar form to the Hamiltonian con- 
straint of the incoherent dust derived in ||| . There is one 
difference, namely that the Hamiltonian constraint 
of the incoherent dust depends on the momenta conju- 
gate to the 3-dimensional coordinate frame variables only 
through the momentum constraint, while in l|63l) Pl ap- 
pears both explicitly and through H 2ND . In spite of 
this, we can still follow the algorithm of as will be- 
come transparent in the following. 

The ADM decomposition of the total action leads to 
the super Hamiltonian and super momentum constraints 



H,. 



--H^ + H\ ND = 
--H? + H 2ND = 



(65) 
(66) 



where the vacuum constraints H^_ and H^f are expressed 
in terms of the canonical variables (M, R; Pj\/, Pr) in Eqs. 
© and lO- By using Eqs. Jj2J}, and (EHJ we can 
eliminate P L , H™ and H 2ND from Eq. JHSJ) to obtain 



= 2\ 



7 (Z' 2 + L' 2 R 2 )P Z + Z'H? 
\ L'AR 



(67) 



Then Pz can be separated from the other variables in 
Eq. (E3 : 



P T z =Pz + h z [M, R, Z, L, P M , P R ] = , 
L'ARh + Z'H? 
Z ~ Z' 2 + L' 2 R 2 ' 



(68) 



From (EHl we know 



Pz 



L'ARh + Z'H? 
Z' 2 + L /2 R 2 



(69) 



Here we used the notation 1121) . 

By using and l|6^|l , the constraint (|5fj|l can be writ- 
ten as 

= H r := + H 2ND 
= + Z'P Z + L'P L 

H?-Z> L '^ h + Z : H / + L>P L , (70) 



Z' 2 + L I2 R 2 



Which gives 



= P L 



-Z'ARh + H^L'R 2 



Z' 2 + L' 2 R 2 
We will denote the constraint l|71|) by H^l. 

H 1L =P L +n L [M, R, Z, L, P M , P R ] = 
-Z'ARh + L'R 2 H? 



(71) 
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Z' 2 + L' 2 R 2 



(72) 



Thus we have obtained a new, more convenient set of 
super-Hamiltonian constraint H\z and supcrmomentum 
constraint H-\l- Both linearized constraints contain ex- 
actly one null dust momentum. The Dirac algebra of the 
old constraints turns into an Abelian algebra of the new 
constraints: 



{H u (r),H ul (r')} = 



(73) 



where H^j = (H^z, H^l) . This feature is similar to the 
case of the one-component ordinary dust |6j , and in fact 
the proof proceeds exactly in the same way. Following 
first we note that the Poisson brackets of the new con- 
straints must vanish, at least weakly (on the constraint 
hypersurface). However, due to the linearity of the con- 
straints (|68|) . I|72l) in the momenta Pz, Pl, the brackets 
do not depend on any of Pz , Pl ■ But then there is no way 
the constraints (|68|l . I)72|l would help in turning into zero 
the Poisson brackets. Therefore they have to strongly 
vanish. 



V. FALLOFF OF THE CANONICAL 
VARIABLES 

A. Falloff conditions for the eternal Schwarzschild 
black hole 

The proof of Kuchaf in Ref. 0] that the mapping 
(A, R, Pa, P r ) — > (M, R, Pm, Pr) is a canonical trans- 
formation in the gravitational sector relies on the check 
that the difference of the Liouvillle forms is an exact form. 
This translates to show that the expression 



B(r) = -RSRln 
y ' 2 



RR! + AP A 



RR 1 - AP A 



(74) 
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vanishes on the boundaries of the domain of integration. 
For the eternal Schwarzschild black hole discussed there, 
the desired behaviour was assured at r — > ±00 by impos- 
ing suitable falloff conditions for the canonical variables, 
based on the treatment of Beig and O'Murchadha |24| . 
The proper falloff of the variables A, R, P\, Pr, Killing 
time T, lapse function N and shift N r , given by Eqs. 
(49)-(55) of Ref. [3 assure that the Kuchaf mapping is 
a canonical transformation. 



as 

H G = -2P A . 2 A r + 0(r 2 ) . (78) 

Thus, the falloff conditions are consistent with the vac- 
uum constraints. They are also preserved in time, as 
noted in pr| . but we will show that only for 

Pa 2 = . (79) 



B. Falloff conditions for r — > in flat space-time 

Hajfcek and Kiefer have studied the evolution of a 
spherically symmetric null dust shell in the space-time 
generated by an other spherically symmetric null dust 
shell The (innermost) region surrounded by the in- 
coming null shell is Minkowski. In order to avoid the 
occurrence of a conical singularity at r = 0, followin g th e 
method developed for cylindrical gravitational waves [25| , 
they have imposed boundary conditions on both the coor- 
dinates and their spatial derivatives at the regular centre. 
Based on these, Bicak and Hajfcek 01 have shown that 
the boundary term l|74[) also vanishes at r — > 0. 

Louko, Whiting and Friedman have discussed the 
Hamiltonian dynamics of a thin (distributional) null-dust 
shell under both sets of boundary conditions: first at the 
two spatial infinities r — > ±00 of the Kruskal-like mani- 
fold and second at r — > and r — ► 00 j2fj. In the latter 
case, the falloff conditions at r — ► for the canonical 
variables, lapse and shift in the flat geometry within the 
null shell are given by their system of Eqs. (7.1): 



A(i,r) 


= A + 


0{r 2 ) , 


R(t,r) 


= R x r- 


KO(r 3 ) , 


P\(t,r) 


= Pa 2 t 


2 +G(r 4 ), 


P R (t,r) 


= P Rl r 


+ 0(r 3 ) , 


N(t,r) 


= N Q + 


0(r 2 ) , 


N r (t,r) 


= N[r - 





(75) 

where A , R\ 1 P\ 2 , Pr^i Af and N[ are functions of 
time. With these falloffs, the expression B (0) vanishes, 
in accordance with the conclusion of Ref. [TJj . 

Given the falloff behaviors (|75|l . all terms in the vac- 
uum gravitational super-Hamiltonian constraint l|18|) are 
O (r 2 ), with two exceptions: R' 2 /2A = i? 2 /2A + O (r 2 ) 
and -A/2 = -A /2 + O (r 2 ). Therefore 



R\ 



2A 



A 2 



The leading term vanishes for 

Ri = A • 



(76) 



(77) 



For this choice, the falloff conditions obey the vacuum 
gravitational super-Hamiltonian constraint. The gravita- 
tional super- momentum constraint i|19|) . in turn behaves 



This can be seen from the following argument The 
time-evolution of the super-Hamiltonian and super- 
momentum constraints are linear combinations of the 
constraints and their covariant derivatives on the leaves: 

Hi = 2HfD r N - 2NKHf 
+ N r D r H? + ND r H Gr , 



Hi = 2HfD r N 



NKPt, 



G 



H?D r N r 



ND r Hf + N r D r H? 



(80) 



Here K = A^Rr 2 (RP R - AP A ) + 2R- 2 P A is the trace 
of the extrinsic curvature of the leaves St., given in Ref. 
[Tlj . From the falloff conditions (|75|l we obtain 



K = A^R^ 2 {R^r, + A P A J + O (r) 



(81) 



Thus the terms proportional to the gravitational con- 
straints, whether they contain K or not, will decay at 
least as O (r) and O (r 2 ) , respectively (provided R\ = A 
was chosen). Problems could aarise only from the terms 
containing derivatives of the constraints. The falloff con- 
ditions (|75|l and the covariant derivatives of the scalar 
and vector densities and H G give at r — > 



N r D r Hf 
ND r H Gr = -2N 



K>2 _ a 2 
N rR L A0 +o{r) 



A 2 



Pa, 

An 



0(r) 



ND r H' 



N r D r Hr 



G 



-N 



R 2 ~A 2 
2A n 



(2r" a 



A n 



0(r) 



= -2NlA P A2 r + O (r 2 ) 



(82) 



Therefore 



jjG = _ N rR L _K _ 2Nq ^ + Q (r) 
A n Aq 



l 



Hi 



An" 1 ) 



0(r) 



(83) 



By chosing the condition l|77|) . the expression for H G 
will decay as O (r). As we exclude the possibility A*o = 
(which would froze time evolution at r = 0), the only 
possibility remaining for a proper decay of H G is to set 
P\ 2 = 0, which completes our proof. 
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Falloff conditions for the radiative atmosphere 
of a star 



Now we study the question, whether the Kuchaf map- 
ping (A, R, P A , P R ) -> (M, R, P M , Pr) of the gravi- 
tational variables remains a canonical transformation in 
the configuration discussed in this paper. In order to 
answer this question, first we remark that the range of 
the cross-streaming null dust metric parameter B is re- 
stricted by R > 0. This determines a lower boundary 
L m in of L, corresponding to R = 0. The quasi-local 
mass function m(L) (for a definition see vanishes 
at some L TO= o(a, B) and takes negative values below, in 
the interval L m i n < L < L m= o- Besides, for L — ► oo the 
solution l|27(l is not asymptotically flat. One can escape 
these unpleasant features by cutting off the space-time 
between certain L± > L m= o and an appropriate high 
value L2 > L\ and matching with appropriate metrics 
across these boundaries (see Fig. ^| . The cross-streaming 
null dust region l|27l) is matched then from the interior 
with the interior Schwarzschild solution, representing a 
static star with mass M\ , whereas from the exterior it is 
bounded by incoming and outgoing Vaidya regions, and 
it touches three exterior Schwarzschild regions in three 
points (these are 2-spheres, if we take into account the 
angles 9, if). Therefore the solution (|27|l is interpreted 
as a thick shell of 2-component radiation, created from 
the intersection of incoming and outgoing thick radiation 
shells. 

The intersection of the last incoming ray with the 
first outgoing ray is the point (2-sphere) where the junc- 
tion to the outermost Schwarzschild region (character- 
ized by mass M2) is done. This region extends towards 
the spatial infinity i . As the fluid region is bounded, 
only the proper falloff at i° of the gravitational variables 
M, R, Pm-, Pr has to hold, as summarized in the first 
subsection of this Section. 

The situation is not so trivial on the other boundary, at 
r — > 0. There, in contrast with the previous treatments 
of Refs. 01 1 13 an d US 1 we do not have vacuum, but 
rather the center of a static star represented by the inte- 
rior Schwarzschild solution, where the falloff conditions 
are not yet known. The linc-clcmcnt representing the 
gravitational field in the interior Schwarzschild solution, 














d 
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n* 


1 

t 








\ 

■ 




7 dusl 
static 




Schwarzschild 
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FIG. 1: (Color online) The static star, incoming and outgo- 
ing radiation zones, crossflowing null dust region and three 
exterior Schwazschild domains. The normal vector n a of the 
foliation leaf Et is different from the fluid 4- velocity u a . The 
foliation is chosen so that time evolution in the static star 
region proceeds in the direction of the fluid world-lines. 



to the ^-velocity u a of the fluid particles.) are given as 



p = const , 

bF 1 ' 2 ~ a/3 



P = P- 



bF 1 / 2 



(87) 



Here k 2 — 8irG and a, b are constants, chosen such that 

P > 0. 

As the canonical treatment of the interior 
Schwarzschild solution has not been developed yet 
(and it is beyond the scope of the present paper), we 
will impose the simplifying condition that the worldlincs 
of the fluid particles of the stellar material are along the 
time evolution vector d/dt. 



ds 2 = - a-bF 1/2 ) dt 2 + F- X dr 2 + r 2 dn 2 



F (r) = 1 



2 2 

k pr 



(84) 
(85) 



is generated by a perfect fluid with energy-momentum 
tensor 



Tab = (p + p) u a u b + p {A) g ab , 



(86) 



where the energy density p and pressure p {with respect 



otu 



—r- = Nn a 



V A- 1 (#- 

or 



where a (t, r) > is a scaling function. From the con- 
dition of normalization of the 4- velocity u a u a = — 1 we 
obtain a 2 = N 2 — N r 2 . This choice of the allowable 
foliations is in accordance with the generic expectation, 
that whenever a reference fluid is present in the system, 
it is advantageaus to introduce the parameter associated 
with the world-lines of the reference fluid as time variable. 
Outside the interior Schwarzschild region, the leaves E t 
are still allowed to be arbitrary space-like hypersurfaces. 
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The energy density and energy current density of the 
fluid with respect to the chosen foliation become 



(I = T ab n a n° =( — ) P+{—) V, 



a h 



j r = T ab n a g. 



NN r 

—^(P + P) 



With the falloff conditions (|75|l at r — > the condi- 
tion B — > will continue to hold, thus the Kuchaf trans- 
formation is canonical. But are these falloff conditions 
consistent with the constraints? In order to responde af- 
firmatively, first we note that for the fluid variables we 
have the following falloff conditions 



36 



P 



Ps + O (r 2 



3 (a - b) 
a = N + 0{r 2 ) , 
p = Ps + 0(r 2 ) , 

> = -JTli^V + 0{r 2 ) . (89) 
3 (a — b) 

These, together with ^Jg = AR 2 sin 9 imply that 



j^j-sta> 



_L 

star 



= 2K 2 Jg l i = 0(r 2 ) , 



(90) 
(91) 



which shows that the total constraints of gravity and fluid 
are obeyed for the chosen falloffs on the boundary, pro- 
vided the condition l|77|l holds. 

The last question to address is whether time evolu- 
tion conserves these falloffs. In order to see this we 
note that both H^ ar and H" tar vanish for the interior 
Schwarzschild solution and for the chosen class of fo- 
liations as 0(r 2 ) and 0(r 3 ), respectively. Thus, the 
falloff of the matter part of the constraints is faster 
than the falloff of the gravitational part, given by Eqs. 
(|S3|l . Fulfilling the conditions J77J and is sufficient 
for the consistency with the constraints in the interior 
Schwarzschild solution. . 

Alternatively, if we do not insist on the interpretation 
of the cross-streaming null dust space-time region as a 
radiation atmosphere of a star, we can let the outgoing 
radiation to emerge from the origin and the incoming 
component to be absorbed by the boundary at r = 0. In 
this case Cauchy surfaces can be chosen in such a way, 
that their boundary at r — » is in a flat space-time, as 
in Fig. 3. of Ref. 14] . In this setup, the expression B 
again vanishes, and the Kuchar transformation is proved 
to be canonical. 



VI. COMPARISON WITH PREVIOUS RESULTS 

In this Section we will establish the connections be- 
tween the sets of canonical variables employed in this 
paper and in Ref. 14]. In order to do this first we il- 
lustrate in Subsec. V.A that an internal time can be 



introduced for a generic spherically symmetric crossflow 
of radiation streams. We start from the variables em- 
ployed in Ref. 14]. In Subsec. V.B we show that the 
connection of those variables with our variables can be 
written up explicitly. 



A. Constraints of null dust crossflow 

Bicak and Hajfcek 0] generalized the canonical for- 
mulation of the one-component null dust, presented in 
Ref. |rjj for a two-component null dust, with the speci- 
fication of spherically symmetry. The gravitational part 
of their action was given by (|17|) . whereas the energy- 
momentum tensor has been chosen as 



T ab = Qap + + jo jb ) j 



with 



1" 



n±<n 



AR 



n a ±A- 1 



(92) 



(93) 



being the four-velocity null vectors of the ingoing and 
outgoing null dust streams. The latter were characterized 
by the canonical coordinates , <f>_ and their conjugate 
momenta n+ , n_. The unit normal to the leaves was 
denoted n a . The canonical action of the system became 

S T [N, N r , A, R, $+, P A , p r , n+, n_] 



dt j dr(P A A + P R R + U + ^+ + n_$_ 

- NHl- N r Hj) , 
with the super-Hamiltonian constraint 



(94) 



HT := H? + Hf H = 



H 



BH 



A 



and super-momentum constraint 

H^ := H r + = , 

H? H = IT $' +n_$'_ . 



(95) 



(96) 



Following the convention of Ref. ^4(, we assume that 
Thus we will use 



n+$^ < o < n_$' 



H BH = 



-n + $' + + n_$'_ 

A 



(97) 



The constraints l|95|) and l|9tj|> can be conveniently com- 
bined as follows 



= 
Similarly 



aht - m 



o = 



AH 



H. 



2$'_ 



-AHf + ff r G 
2<P' 



AH? + H? 



2$'_ 



(98) 



(99) 
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The new constraints (I98|l and (|99(l are analogous to the 
previously introduced constraints Hti8[) and l|72l) in con- 
taining the canonical momenta only linearly. 

There is nothing to prevent us in introducing square- 
root type new constraints by properly transforming the 
constraints l|98|) . I|99|) as well. The product of the null 
dust momenta is 



n+ir 



- (AH 



BH\ 



BH\ 



4<j ) /2<j ) /2 



(100) 



where we have employed the constraints ()98|l and Ij99(l 
in the first equality and the constraints Q95JI and l|96|) in 
the second equality. Now we can express either H± H or 
H^ H as a square root: 



H 



BH 



n+n_ + g rr HBH H BH t 

A 



BH\ 



h? h = ^/$^_$'_n + n_ + (ah 

These formulae are analogous to the result I|tj3fl . 
The simple linear transformation 

T=i= + , 

.= i= 

p T = -^(n + + n_) , 
P. = ^(n + -n_) 



(101) 



(102) 



B. Static case 

In this subsection we establish the connection be- 
tween the canonical variables (Q±, H±) employed in Ref. 
[l^ | and our canonical coordinates (Z. L) and momenta 
(Pz, Pl)- hi order to do this, first we introduce the 
"tortoise- type" radial coordinate R* defined as dR* = 
R (L) dL. Next we define null coordinates X± — Z ± R* 
(see Ip). 

As in the static scenario the metric is uniquely given 
by Eq. I|27l) . we would like to identify the energy- 
momentum tensors (|92H and (|25|) . which yields — 
K 2 u a /R 2 W and l a _ = n 2 v a /R 2 W, with a proportional- 
ity constant k 2 . According to Ref. [14| the null forms 
l± = A 2 (|II±| / yfg) d&± (A 2 a possible second propor- 
tionality constant). By employing Eq. H29fl . we conclude 
that 



( K 
A 



dX± 



(105) 



with R regarded as a function of X± . 

Equivalently, the derived canonical coordinates (T, a) 
are related in a simple way to (Z, L): 



^{2k\ dZ 
P a dT + P T da = -V? ( — ' 



R 2 



P T dT + P„da = V9 (f) f 



(106) 



with R representing the function l|28[) of L. 

The transformations (|105|l and (|106|1 establish the re- 
lation between our results and the results derived in Ref. 

m. 



is a canonical transformation as can be checked by calcu- 
lating the Poisson brackets of the new canonical variables 
T, a, P T , P a ■ The sum of the constraints (JHHJl and (|9~9"j) . 
divided by V2, in the new canonical chart becomes 



= P T 



a' AH? + T'H? 



(103) 



Similarly the difference of and , divided by 
reads 



—T' AH? - a'Hf 
n — p -I ± L. 



(104) 



The canonical momenta of the cross-flowing null dust are 
then completely separated from the rest of the variables 
in the new constraints (|103fl and (|1U4|> . 

If the null vectors l± a oc $±, a are both future- 
oriented, then T is time-like coordinate. Therefore the 
quantization of the constraint i|103|) will give a functional 
Schrodinger equation. 



VII. CONCLUDING REMARKS 

The geometrical optics approximation of radiation 
fields is represented by null dust. This approximation 
is a very good one whenever the wavelength of the radi- 
ation is short as compared to the typical local curvature 
scale of the space-time. The crossflow of two such ra- 
diation streams describes the interesting situation of a 
two-component radiation atmosphere of a star. The as- 
sumptions of spherical symmetry and staticity lead to 
the space-time (|27|l . The two null dust components in- 
teract only gravitationally (through the curvature of the 
space-time they jointly produce) and formally they are 
equivalent to an anisotropic fluid. 

A previous canonical treatment of such a system |14| , 
besides its many achievements, still suffers from the lack 
of an internal time (a difficulty first encountered in the 
case of one-component null dust) . The existence of such a 
preferred time function would simplify the quantization 
of the gravitational field in question. The absence of 
the internal time from the formalism of Ref. Il4l is not 
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a major inconvenience for the analysis presented there, 
where the one-component null dust limit (Vaidya space- 
time) is discussed in detail. 

If one does not aim to have this limit in the formal- 
ism, the situation is different. We have shown on the 
example of the static crossflow of radiation streams how 
to construct an internal time for the two-component sys- 
tem. By suitable canonical transformations we have in- 
troduced the time function Z as canonical coordinate 
and we have constructed the new super-Hamiltonian and 
super-momentum constraints, Eqs. (|68[) . (|72[) . which 
have strongly vanishing Poisson brackets. With this, we 
have turned the Dirac algebra of the original constraints 
into an Abelian algebra. 

The new constraints contain the momenta conjugate 
to the crossflowing null dust variables linearly. This 
convenient feature can be further exploited in the pro- 



cess of quantization, which will turn the new super- 
Hamiltonian constraint into a functional Schrodinger 
equation. The latter has the obvious advantage over the 
Wheeler-deWitt equation obtained by the quantization of 
the original super-Hamiltonian constraint, that its space 
of solutions is linear. Further properties of the resulting 
functional Schrodinger equation are under investigation 
and we propose to discuss this topic in detail elsewhere. 
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